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Abstract 

Point polarizable molecules at fixed spatial positions have solvable electrostatic properties in classical 
approximation, the most familiar being the Clausius-Mossotti (CM) formula. This paper generalizes the 
model and imagines various applications to nanosystems. The behavior is worked out for a sequence of 
octahedral fragments of simple cubic crystals, and the crossover to the bulk CM law is found. Some 
relations to fixed moment systems are discussed and exploited. The one-dimensional dipole stack is 
introduced as an important model system. The energy of interaction of parallel stacks is worked out, 
and clarifies the diverse behavior found in different crystal structures. It also suggests patterns of self- 
organization which polar molecules in solution might adopt. A sum rule on the stack interaction is 
found and tested. Stability of polarized states under thermal fluctuations is discussed, using the one- 
dimensional domain wall as an example. Possible structures for polar hard ellipsoids are considered. An 
idea is formulated for enhancing polarity of nanosystems by intentionally adding metallic coatings. 

1 Introduction 

Electrical polarity (neutral charge distributions with positive and negative charge centers separated in 
space), and in particular dipole moments, are ubiquitous in asymmetric molecules, and therefore no 
doubt very common in nanocrystals [1]. Systematic studies of polarity of nanocrystals are rare. Two 
primitive questions arise. (1) When can a nanosystem develop a spontaneous polarity, analogous to a 
bulk ferroelectric, i.e. reversible in a strong applied field, and how can this polarity be controlled [2]? A 
particularly fascinating example is the recent discovery [3] of spontaneous polarization in small clusters 
of Nb atoms at low temperature T. (2) What "permanent" moments occur, how stable are they in time, 
and how much fluctuation occurs within an ensemble of similarly prepared nanocrystals? 

The venerable "Clausius-Mossotti model" (CM model, or CMM) [4, 5, 6], provides a simple solvable 
picture of spontaneous polarity. Most of this paper consists of imagining possible applications of the 
CMM and working out some consequences in nanosystems. I define the CMM to mean a system of fixed 
points fi, each of which has a polarizability oq. When an electric field i^tot appears at point i, it induces 
an electric dipole moment fa = o^i^ tot- Each such moment creates a dipole field. The sum of all such 
dipole fields is the induced field 
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and the total field is the sum of these plus the external field -F^tot = -Fi,cxt + Fi,ind- This problem requires 
self-consistency to find the answer for 

The limitation to induced moments and dipole-dipole interactions inhibits application of the CMM to 
real systems. However, by focusing on the dipole-dipole interaction, this model clarifies a difficult aspect 
of real polar systems, including systems with fixed moments. 

There is a large chemical literature developing and using model interactions between molecules, in 
which some or all of the quantum electron behavior can be represented approximately as classical charges 



and dipoles which interact electrostatically [7, 8, 9, 10, 11] Ideas of this kind are particularly appropriate 
for interacting organic molecules (e.g. molecular crystals [12, 13, 14]). Interactions alter their electronic 
properties from values found in isolated molecules; the effect is large, but not so large as to destroy the 
identity of the molecule as happens in covalent, ionic, or metallic interactions. The alterations arise from 
the polarization of surrounding molecules and can be modelled by classical electrostatics. The present 
paper is motivated by more primitive questions; the methods developed here can be found in related but 
more sophisticated contexts elsewhere. 

2 The "Generalized Clausius-Mossotti Law" (CML) 

The textbook version of the CML considers an infinite periodic medium with polarizable molecules sitting 
on equivalent sites of cubic symmetry. The result is that the total field F at any molecular site is enhanced 
from the value ^ oxt of the homogeneous applied field by an enhancement factor / = (1 — 47ma/3) _1 , 
where n is the number density of the molecules. Equivalcntly, a is enhanced to the value fa by the fields 
of the other dipoles. A general derivation from microscopic response theory was given by Maksimov and 
Mazin [15]. 

A good example is crystalline Cqo which has weakly-coupled molecules of known [16, 17] polarizability 
a=85A 3 sitting on a face-centered cubic (fee) lattice of density n—0. 001406 moleculcs/A 3 . The enhance- 
ment factor is / = 2.00 and the dielectric constant [18] is = 1 + Aimfa = 4.0. These independently 
measured values of the the vapor-phase property [16] a, and bulk-solid property [18] e^, provide an 
excellent confirmation of the CML. Metallic screening by intermolecular charge transfer, while needed in 
a complete theory, is small enough in Ceo that the Clausius-Mossotti approach works well. 

The microscopic calculation of Pederson and Quong [17] shows that, for Ceo, the theoretical free- 
molecule polarizability (83. 5 A 3 ) is strongly affected by intramolecular screening. Even within a single 
molecule the local field is reduced by a factor of three by intramolecular charge rearrangement. But when 
assembled into the fee solid, the details of local charge are successfully ignored by the CMM which places 
all the polarity at a single point of cubic symmetry. There is no guarantee that this would be equally 
successful for less symmetric molecules. Schemes for representing the quantum electrostatics of molecules 
by interacting "submolecules" have been developed [12, 13, 14]. 

Note that the CM enhancement factor / diverges to infinity when the polarizability increases to 
a m ax = 3/47rn, or alternately, the density increases to n max = 3/47ra. There is no known molecular 
crystal where this limit is reached, but Cgo comes within a factor of 2. If the limit were reached, the 
material would enter a phase of spontaneous polarization. Actual solids which do polarize spontaneously 
[19, 20] are more complicated, but the CMM provides a useful simple system with the same property. 

To derive a more general result, note that the dipole moments {&} or their Cartesian components 
{A*ia}> obey a minimum principle. That is, they minimize the energy expression 

£di P ({&» = e (I; " A ■ ( H ind + ^'° xt) ) • (2) 

The first term expresses the fact that it costs less to polarize a molecule when its polarizability is large. 
The factor of 1/2 in the second term avoids double counting the dipole-dipolc interaction. The energy 
formula is a quadratic form, amenable to methods of linear algebra. It is helpful to use a vector space 
notation, 

£di P ({&}) = ;r < Ml-i - T\n > - < Ml^cxt >, (3) 
2 a 

where the SA^-dimensional column vectors have been introduced: 
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The 3iV x 3-ZV-dimensional matrix T is the dipole-dipole interaction, defined by Eq.(l), i.e. [find > 
T\fj, >, and has elements 

iRij a Rijfj — SapRf 

~r1. 
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Ef^- r ) (6) 



where Rij a is the a Cartesian component of the vector Rij = Ri — Rj. The solution for the dipoles jli 
which minimize the energy is 
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This is the generalized Clausius-Mossotti law (CML), valid for arbitrary collections of polarizable neutral 
molecules [21]. For simplicity it is assumed that all molecules have the same polarizability. The general- 
ization to multiple species is straightforward. The structure n « faF ext is evident. In sec. 4 it will be 
shown how the usual textbook CML follows from this formula in the appropriate limit. 

It is interesting and useful to consider the eigenvalues 7 and eigenvectors I7 > of the dipole-dipole 
interaction V: 

r = 5>l7><7l- (?) 

7 

In terms of these eigenstates, the generalized CML for the induced dipoles is 



In the usual case, the factors (1/a — 7) _1 are all positive finite numbers. The largest eigenvalue 7 is less 
than the resistance 1/a to polarization. Because the interaction Eq.(5) scales as 1/R 3 , the eigenvalues 
7 scale linearly with density n. If one imagines increasing the density of the system, leaving everything 
else fixed, eventually the largest eigenvalue 7 max exceeds 1/a and the response is divergent. This signifies 
that a spontaneous polarization now occurs, with a pattern of polarization given by the corresponding 
eigenvector |7 max >. 



3 Octahedral Fragments 

How does the CMM describes the evolution from nanocrystal to bulk ferroelectric? Consider the sequence 
of octahedral fragments of a simple cubic lattice. The simplest fragment with a central molecule and 
six octahedral neighbors is shown in Fig. 1. Larger fragments have 7V=25, 63, ... molecules and retain 
the full cubic point symmetry. The largest eigenvalues of the dipole-dipole matrix have been computed 
numerically for these cases and are shown plotted versus a convenient reciprocal power of iV in Fig. 1. 
The N = 00 bulk limit is nearly achieved for 200-molecule clusters. For the infinite simple cubic crystal, 
the CML tells us that a ferroelectric instability occurs when 1/a diminishes to 4tt/3. For fragments of 
the crystal, the critical value of 1/a should surely be less than 47r/3, increasing to this value as N — > 00. 
Fig. 1 shows that there is a monotonic increase with N but that already by N=25 the largest eigenvalue 
exceeds the bulk CML value 47r/3. 

The reason for this unexpected result is that the sc structure prefers antiferro rather than ferro 
ordering when only dipole-dipole interactions are present. This was discovered by Luttinger and Tisza 
[22]. The next section elaborates, and Sec. 7 gives a simple way to calculate the preferred dipole ordering 
on a lattice. To confirm the antiferro preference, Fig. 2 shows how the polarization directions evolve as 
the octahedral fragment enlarges toward the bulk limit. 

A recent report by Fu and Bellaiche [23] considers dipoles on larger cubic fragments of a simple cubic 
lattice. They find a surprising spiral pattern to be stable, with no net moment. They are solving a 
somewhat different model. 



3 



6 




Figure 1: Largest eigenvalue of dipole-dipole matrix for octahedral fragments of a simple cubic lattice. 
Units are used where the density n=l. The smallest octahedral fragment with 7 molecules is shown. The 
horizontal line is the CM eigenvalue 47r/3. 



4 Bulk Clausius-Mossotti Law 

Infinite periodic crystals have translational symmetry that permits computation of the eigenvalues and 
eigenvectors of the dipole-dipole matrix T. Specifically, the Floquet theorem says that eigenstates I7 > 
can be chosen to have the Bloch form |fcA > 

<ia\kX>=e i ^ Ri e a {k\)/VN, (9) 

where (for a one-atom cell) the eigenvector e(fcA) is a 3-d unit vector, labeled by wavenumber k and 
branch index A. This eigenvector corresponds to eigenvalue "f(kX) and satisfies 

T aP (k)ep{k\) = 7 (fcA)e Q (fcA), (10) 

where the 3x3 dipole-dipole matrix T is now in Fourier space, 

For the sc structure, the eigenvalues ^(kX) are shown in Fig. 3. The numerical calculation used the 
Ewald method to converge the sums [24]. The sum of the three eigenvalues (trace of the matrix Ta^^k)) 
is 0, independent of k. This is because V 2 (l/r) is everywhere except at r = and the Ri = Rj term 
is omitted from Eq.(ll). The same dispersion curve shown in Fig. 3 has a different interpretation in 
phonon physics. By a change the sign and an additive constant 47r/3, these become the curves of squared 
frequencies of vibration of the Coulomb lattice of unit point charges. These are the "bare phonons" 
encountered in older approaches to the theory of phonons in metals [25]. The fact that the eigenvalue 
at the r point (k — 0) is not the extreme eigenvalue in sc structure corresponds to a linear instability 
(imaginary eigenfrequencies) of the Coulomb lattice in sc structure. Conversely, the fact that the k = 
eigenvalue 7 is maximal in fee and bec structures corresponds to linear stability of the Coulomb lattice 
in these structures. 

The Clausius-Mossotti formula (/i^ = faF cxt with 1// = 1 — Airna/S) can be obtained from these 
results. The external field is homogeneous, so only eigenvectors | 7 > corresponding to k = can couple 
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Figure 2: Eigenvector properties corresponding to the largest eigenvalue of the dipolc-dipolc matrix for 
octahedral fragments of a simple cubic lattice. The "polar" sites (immediately above and below the 
central ^-polarized site with polarization /i = 1) have their polarization parallel and evolving toward 1 
at large N. The "equatorial" sites (four sites surrounding the central site in the xy plane) have their 
polarization antiparallel and evolving toward -1. 



in Eq.(8). Then from Eq.(8) we get 

3 /i _ \-i _ _ 

rt = Xj(--^( fc = ' A )J (c(fc = 0,A)-F ext )e(fc = 0,A). (12) 

Thus it is needed that all three eigenvalues j(k = 0, A) should be 47m/3. Two eigenvalues indeed agree 
with this, but consistent with summing to zero, the third eigenvalue is — 87m/3. The corresponding 
eigenvector e(k = 0, A) is "longitudinal," that is, it pointed along k before the k — > limit was taken. 
Longitudinal fields suffer depolarization which transverse fields avoid. A longitudinal field is created by 
a parallel plate capacitor applied to a block of dielectric. The "external" field is related to the classical 
macroscopic field F by F ext = F + AirP and the polarization density P is njl. Thus for a longitudinal 
field, Eq.(12) is the same as 

M=(^ + ^p) \p + ^n{l). (13) 
This is equivalent to /x = faF, which is the correct Clausius-Mossotti law for a macroscopic sample. 



5 Fixed rather than Induced Moments 

The CMM refers to a system of polarizable molecules whose induced dipoles ft minimize the energy Eq.(3). 
If spontaneous moments are just barely stable in no applied field, the solution is [i ia oc< 2a|7 max >, where 
l7max > is the eigenvector of T with largest eigenvalue. The proportionality constant is determined by 
higher order terms. 

Another model, more in touch with chemical reality, has molecules with moments fli fixed in magnitude 
= Mo) but free to choose their spatial orientation. The ground-state moment configuration then 
minimizes Eq.(3), except with the resistance 1/a to polarization set to zero, and with the N constraints 
= /io- In general, this is a difficult algebraic problem. However, suppose the eigenvector |7 max > 
has the property that = where /ii a cx< ia\-f max >. For example, an infinite system with perfect 
ferroelectricity or simple antiferroelectricity has this property (each moment has the same magnitude.) 
Then this also must provide the ground state of the fixed moment problem. As a magnetic example 
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Figure 3: Wavevector dependence of the eigenvalues j(kX) of the dipole-dipole operator computed for the 
simple cubic structure, in units where n — l/a 3 = 1. Note that the largest eigenvalue is at K = (it, it, 0). 
At T = (0,0,0), the eigenvalues are 4tt/3 (the Clausius-Mossotti value, doubly degenerate) and — 8ir/3 
(singly degenerate). 

(magnetic and electric fixed moments obey the same dipole-dipole laws) , White et al. [26] studied crystals 
of hydrated rare-earth phosphomolybdates. The magnetic moments of the rare-earth ions (located on 
a diamond-structure sublattice) are fixed by Hund's rules but sufficiently separated from each other 
that interatomic exchange is ignorable. Then the magnetic order (seen below an ordering temperatures 
< 50 mK) is determined by dipole-dipole interactions, and has been measured. The magnetic order is 
antiferromagnetic, much like the result of Sec. 3 for sc structure, except the order is more interesting in 
diamond structure. As explanation, the authors computed the largest antiferromagnetic eigenvalue of the 
matrix T and showed that for diamond structure (and also for sc but not for bec or fee) this eigenvalue 
exceeds the ferromagnetic (CM) eigenvalue 47rn/3. Their calculation is the N — > oo limit of Fig. (1). 

In general, as illustrated in Fig. 2, but unlike the bulk examples just discussed, the vector |7 max > 
with largest eigenvalue of T does not have all moments [ii a oc< ia\-f ma ^ > equal in magnitude, and 
thus does not solve the problem of finding the most favorable alignment of fixed dipolc moments. In 
what sense is the solution of the fixed dipole problem {^} "close" to the unconstrained solution given 
by |7max >? It is required to minimize the dipole energy — < /i|r|/i > /2 plus constraint terms Xi/sf/2. 
The solution depends on the Lagrange multipliers (LM) Xi which are then varied until the constaints 
/if = const arc satisfied. There are N moments pi, but only N — 1 constraints, because the constant 
value of /if is arbitrary since the problem is linear. The minimization can be formulated as a "generalized 
Hermitean eigenvalue problem." The N LM's define a vector A in the A^-dimensional LM space. Write 
this as A = Ac where the components Ci of c are normalized by c i — N. Also write /if =< /-i\Pi\fi > 
where Pj is a projection operator onto the 3- vector space of the i-th moment fli. Then the formal solution 
of the minimization problem is 

I> > -XP(c)\n >= (14) 

where the Hermitean operator P(c) = C %P% is the sum of projection operators. The normalized vector 
c defines a "ray" in LM space. There is one special ray Co, where Cj = 1 for all i, so that P(co) is the unit 
operator. For this special case, Eq.(14) is the ordinary eigenvalue equation for T. For other rays, as long 
as no component Cj is zero, P is invertible and Eq.(14) is a generalized Hermitean eigenvalue equation. 
The N eigenvalues X K depend on the direction c of the ray, and define N surfaces A = A K (c)c in LM 
space. Each surface k intersects the special ray at distance j K , where j K is an ordinary eigenvalue of T. 

All that has been accomplished so far is that the eigenvalue problem for T has been generalized to an 
infinite family of eigenvalue problems, one for each ray c in LM space. Now we have to find the optimal 
solution |/i > which minimizes dipole energy with fixed components [n — fi. The N — 1 constraints 
/if = const will be satisfied along 1-dimensional lines in LM space. These lines will puncture each surface 
k at isolated points. The optimal solution occurs at such a point on the surface of maximum eigenvalues. 
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That is, it lies on the surface which passes through 7 max c , i.e., the surface which evolves from the solution 
of the unconstrained problem. This suggests an algorithm. First find the optimal solution |7 max > of 
the unconstrained problem. Then find the nearby eigenvalue A max sa j max and corresponding eigenvector 
|A ma x > of the generalized problem for rays c near the special ray c . Wander stepwise through nearby 
rays guided by the demand that the components jli = Pi|A max > of the maximal eigenvector should evolve 
toward equal lengths. This should provide an efficient iterative solution of the fixed moment problem. 



6 1-Dimensional Dipole Stack 



Polar molecules in solution [27] may tend to self-organize into structures which optimize dipole-dipole 
interactions. The "Stockmayer fluid" [28] is a simplified model for such a system, which keeps hard 
core repulsion and long-range dipole-dipole interactions for spherical molecules with point dipoles. Such 
models may show linear arrangements of dipoles [29] . Such arrangements have also been directly visualized 
for magnetic nanoparticlcs in solution [30], and even for Au nanoparticles [31] which are presumed 
to acquire their polarity from interactions with species on their surface or in solution. This section 
considers ordered one-dimensional arrays, and the subsequent two sections considers the interactions 
between different ordered linear arrays. 

Consider a one-dimensional stack of dipoles, with repeat distance 6. For simplicity, think of them as 
permanent rather than induced moments, with identical magnitudes. Let each moment jl point along the 
stack axis, chosen as the z axis. The total interaction energy is 
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Thus adding a molecule to the stack gains energy 2£(3)/i 2 /6 3 where £ is the Riemann ^-function and £(3) 
has the numerical value 1.202057. . .. Breaking a stack into two smaller stacks creates two new surface 
points, each costing energy (/z 2 /6 3 )£(2) where £(2) = 7r 2 /6 = 1.6449 . . .. 

If the molecules are instead polarizable (with zero fixed moment) then the stack will spontaneously 
polarize when the energy gain exceeds the cost fi 2 /2a, that is, when 1/a is less than 4£(3)/6 3 . One 
should ask if this situation is thermodynamically stable. For example, how much does it cost to rotate 
the direction of jt away from the stack axis? How much does it cost to create a defect, for example, a 
stacking fault where the direction of jl changes at the fault? 

For the octahedral fragments considered in Sec. 3, cubic point symmetry has the consequence that 
until terms of higher order are considered, there is no preference for any particular direction of jl. There 
is no anisotropy energy, so the ordered moment lacks a permanent direction in space. A collection of 
such nanoparticles would be called a "superparaelectric," the analog of a "superparamagnet." For the 
1-d stack, this is no longer the case. Cubic symmetry is maximally broken, and it costs 3£(3) sin 2 (6*)/Lt 2 /63 
per molecule to rotate all dipoles. The anisotropy energy is very large. 

On the other hand, it is a well-known property of ordered states in one dimension that domain wall 
defects destroy long-range order. We can estimate the cost of a domain wall as follows. An absolutely 
abrupt boundary, where induced moments have their full positive value on one side of the boundary 
between two molecules, and their full negative value on the other side, costs energy 4£(2)^ 2 /6 3 in lost 
dipole-dipole attraction. A more gradual case is an abrupt boundary on the site of a molecule. This 
loses the same amount of dipole-dipole attraction but there is a savings p? /2a since one molecule is 
unpolarizcd, making this domain wall structure more favorable. More gradual boundaries become lower 
in energy when the polarizability is diminished. 

Domain-wall defects will always be thermally excited except in shorter chains at lower temperature. 
The number of thermally excited domain walls is Np where p = cxp(— A/fc^T) and A is the energy to 
create a domain wall. Suppose molecules in solution interact by dipolar forces. It is argued in Sec. 8 
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that for oblate ellipsoidal molecules with dipoles along the short axis, chains are favored and only weakly 
attract neighboring chains. Then the statistical distribution of chain lengths n may be given by the 
Flory-Schulz [32] distribution P n = p(l — p)™ -1 . The mean chain length would be 1/p but the most 
probable chain length would be 1. The derivation of this distribution assumes that A and thus p do not 
depend on length, whereas actually A becomes smaller for short chains, which will enhance P n at small 
n. The Flory-Schulz distribution is usually written in terms of a variable p' = 1 — p which measures the 
"completion" of the "reaction." For prolate ellipsoids, there will be greater organization of chains into 
bundles. Domain walls then become correspondingly larger in area and harder to excite thermally, and 
the system evolves into a ferroelectric. 



7 Interacting Parallel Stacks 



There is an interesting and important interaction between two parallel stacks of fixed dipoles. Consider 
the case of infinite stacks (vertical separation of moments is 6), where each moment is vertical (along the 
stack) and the two stacks are identical and parallel. Then it is obvious that the interaction is a periodic 
function of the vertical offset y with period b. It turns out that the interaction is strictly oscillatory 
with no constant component (in powers of exp(i2iry/b)) and decays exponentially (as exp(— 2-Kx/b)) with 
horizontal separation x. 

The stack-stack interaction (per molecule) is 



U K 
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53 Z_, [ x 2 + ( n6 _ y )2]3/2 
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Hnb-yf 



+ (nb - y) 1 



(15) 



This is also the interaction energy between a single infinite stack and a single z-oriented dipole at (x,y). 
At large separations x we can use the continuum approximation of a constant linear dipole density ji/b 
along the stack. Then the stack is electrostatically equivalent to charges ±/i/6 on the ends of the stacks 
which arc infinitely removed from the observation point. Therefore at large separations the electric field 
and the interaction are zero. Any non-zero remainder (contained in Eq.(15) derives from the grainincss 
of the dipole density. 

To derive an asymptotic formula for the interaction, use the Poisson sum formula 



00 ^ 

n — — rvi ^ 



COth(n z)f(z). 



(16) 



The contour in the complex z-plane surrounds the imaginary axis where coth(nz) has poles at z = in. 
Rewrite the summand of Eq.(15) as —d 2 /dy 2 [x 2 + (nb — y) 2 ]^ 1 ^ 2 - After bending the contour around 
suitable branch cuts, the dimcnsionlcss interaction is 
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(17) 



From this exact representation, extract the leading term in the asymptotic expansion in powers exp(2irnx/b) 
and (x/b)-" 1 - 1 / 2 , 

= -27rx/b 



(f^^cos^/fc)- 



(18) 



An understanding of formulas like this dates back to Madelung [33] . Because of y-periodicity, the electro- 
static potential <j>(x, y) can be expanded as ^2 Ce.fi(x) exp(2iri£y/b). Madelung points out that Poisson's 
equation demands that fi(x) be the Hankel function Ko(2ir£x/b). The absence of the constant (I = 0) 
term from the sum follows from the charge neutrality of the stack. The asymptotic behavior of Ko(t) is 
exp(— t)^iT/2t. In Figs. 4 and 5, the exact and the asymptotic version of the interaction are compared. 
They agree to within ±15% for (x/b) = 0.5 and to within ±2% at (x/b) — 1. A plot of the electric field 
near the stack is shown in a separate paper [35] 

The oscillation in u ss as a function ofy/b shown in Fig. 5 has interesting consequences. When adjacent 
stacks have zero offset (as happens in the sc crystal structure) , their moments jl prefer to be antiparallel, 
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horizontal separation x/b 



Figure 4: Stack-stack interaction between two identical stacks of fixed dipoles \i separated along the stack 
by distance b. The stacks are separated horizontally by distance x shown horizontally in units of b. The 
vertical offset y is zero. The interaction is in units ^ 2 /6 3 . The solid line is computed numerically from 
Eq.(15) and the circles are the first term, Eq.(18) in the asymptotic expansion. 

as was already seen for the octahedral fragments in Sec. 3. When adjacent stacks have a full offset 
y/b = 0.5, their moments prefer to be parallel, which is known to be the preference in the bec and fee 
crystal structures. Triangular packing of adjacent stacks has "frustrated" interactions. If the offset is 
zero, then some nearest neighbor interactions are necessarily repulsive no matter what orientations of fx 
are chosen. This problem is discussed further later in this section and in the next. 

For cubic structures, an interesting sum-rule on stack-stack interactions is implied by the CM formula. 
This formula is equivalent to the statement that the dipole-dipole interaction energy for a cubic arrange- 
ment of dipoles, all oriented in the same direction, is — 7cmM 2 /2 per molecule, and the CM eigenvalue 
7cm is 47m/3. Crystal structures can be viewed (in multiple ways) as parallel lines of stacked molecules. 
If the dipoles line up along the stacks, the total dipole-dipole energy can be computed as the energy per 
molecule in a stack, plus the energy per molecule of stack-stack interactions. Set these two versions of 
the total energy equal to each other, divide out /i 2 , and multiply by the volume of a cube whose edge 
length is the periodicity b on the stack. Define v = nb 3 . Then the dimcnsionlcss energy equation is the 
sum rule 

^ stacks 

Y = - 2 C(3) + 2 "ss (y, j) (cubic structures) (19) 

Because the stack-stack interaction decays rapidly with separation, the dipole-dipole energy comes 
from the single stack energy — 2£(3) = —2.4041 plus near neighbor stack interactions (the last term in 
Eq. 19, with a factor of 1/2 to avoid double counting.) This is demonstrated to be true in tables 1- 
3. For each crystal structure, there is more than one simple axis which can be chosen as a stack axis. 
In bee structure (lll)-oriented stacks are a favorable choice relative to (100) orientation because the 
spacing of molecules is closer (b = \/3a/2 instead of a); (110) stacks are best for both diamond and fee 
structures, while (100) stacks are best in sc structure. The sum rule applies for any consistent choice of 
stacks. First-neighbor interactions converge the sum rule to better than 1% for sc (100) and bec (111) 
stacks, and fee requires only second-neighbor interactions in either (100) or (110) orientations. It is easily 
seen how the sc structure lowers its dipole energy by switching orientation of first-neighbor (100) stacks. 
For diamond, although (110) stacks give only a small improvement in convergence, they allow a better 
understanding of the antiferroelectric solution which gives the largest eigenvalue 7. Each (110) stack has 
one very close neighbor with offset y/b = 1/2 which makes a strong ferroelectric bond, and two quite 
close next neighbors with zero offset, which therefore prefer to be antiferroelectrically oriented. 
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0.1 0.2 0.3 0.4 0.5 
vertical offset y/b 

Figure 5: Stack-stack interaction between two identical stacks of fixed dipoles [i separated along the stack 
by distance b. The stacks are separated horizontally by distance x/b — 0.5 (solid line and filled circles) 
or by distance x/b — 1.0 (dashed line and open circles. The vertical offset y/b is shown horizontally. The 
interaction is in units /j 2 /b 3 . Note that the x = b interaction is shown after multiplication by 20. 

Solutions of polar molecules of ellipsoidal shape may self-organize into arrays of stacks. For cylindri- 
cally shaped molecules, the stacks may organize in a 2-d triangular lattice. Therefore it is interesting to 
look at the cubic structures from this point of view. Diamond is the only cubic structure that does not 
work, because (111) stacks are dimerized (alternate spacings of -\/3a/4 and 3-^0/4. occur). The other 
cubic structures are triangular arrays of (111) stacks. The "frustration" problem is solved by alternating 
offsets of ±1/3 for the six first-neighbor stacks surrounding a central stack. At this offset, ferroelec- 
tric alignment is preferred by all first-neighbor pairs. Consider the rhombohcdral translation vectors 
Ri = R{6i) =rcos6iX + r sin^y + (b/3)z, where 0* = 0°, 120°, and 240°. The sum of these three primi- 
tive translations is bz, the period of the stack, and the distance r is the transverse separation of stacks. 
The rhombohedral angle a is defined by R 2 cos a = Ri-Rj which gives cos a = [(6/3) 2 — r 2 /2]/[(6/3) 2 +r 2 ]. 
The three special values of cos a arc -1/3 (bee), (sc), and 1/2 (fee), which correspond to ratios r/b of 
spacing to period of r\ = \/8/9 (bec), ri/2 (sc), and ri/4 (fee). For bee, the period and separation are 
nearly the same, so stacks interact only weakly. The sc structure is intermediate, while fee structure 
has a small transverse separation relative to periodicity, leading to very strong stack-stack interactions. 
For both fee and bee structures, the ferroelectric alignment is optimal. However, for the intermediate 
case of sc structure, a lower energy solution tilts the axis of polarization of molecules along (111) stacks 
alternately to ±(100). Structures like sc and fee should not necessarily be considered as arrays of (111) 
stacks even though it may be mathematically interesting to do so. 
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structure 


V 


-™ + 2C(3) 


Xj 






\MjU ss (xj,yj) 


sc 


1 


0.3097 


1 





4 


0.2910 


bcc 


2 


-1.7847 


i 

1 


1 

2 


4 


-2.0822 











4 


0.2910 














T =-1.7912 


fee 


4 


-5.9735 


1 

2 

1 
V2 


1 
2 



4 
4 


-8.2484 
2.2272 














>J = -6.0212 


diamond 


8 


-14.3510 


1 

1 


1 

4 
1 
2 


4 
4 


-8.3287 
-8.2484 








I 2 
V2 

VTo 

4 





4 


2.2272 








1 
4 


8 


-0.0481 








1 





4 


0.2910 








3 

2 


1 

4 
1 
2 


4 
8 


-0.0007 
-0.2606 














V = -14.3683 



Table. 1 Sum rule satisfaction for cubic crystals where the stacks have period b — a along a cube axis, 
and dipoles are oriented along the same axis, v is the number of molecules in a cube of edge b. Mj is 
the multiplicity of the j-th set of neighboring stacks. 



structure 


V 


=tF + 2 C(3) 


Xj 


Vi 






fee 




-0.5578 


V3 
2 


1 

2 


4 


-0.7110 








1 





2 


0.1455 














>J = -0.5655 


diamond 




-3.5196 


1 

2^2 


1 

2 


1 


-5.0847 








V6 
4 





2 


2.2235 








2 


1 

2 


4 


-0.7710 








1 





2 


0.1455 








3^2 
4 


1 

2 


3 


-0.1436 








\/22 
4 





2 


0.0454 














J2 = -3.5249 



Table. 2 Sum rule satisfaction for cubic crystals where the stacks have period b = a/y/2 along a cube 
face diagonal, v is the number of molecules in a cube of edge b. 



structure 


V 


^P + 2((3) 


Xj 


Vi 


Mj 


±MjU as {xj,yj) 


bcc 


3^ 
4 


-0.3166 


v/879 


1 
3 


6 


-0.3224 


sc 


3^3 


-8.4787 


V 2 / 9 


1 

3 


6 


-9.6448 








V^/9 





6 


1.5423 










1 
3 


6 


-0.3224 








v/14/9 


1 

3 


12 


-0.0826 














J2 = -8.5075 



Table. 3 Sum rule satisfaction for cubic crystals where the stacks have period b = V3a/2 (bcc) and 
b = \/3a (sc) along a cube body diagonal, v is the number of molecules in a cube of edge b. The sum 
rule is satisfied to similar accuracy for fee structure with (111) stacks only when stacks out to 12th 

neighbors are summed. 



11 



8 Close Packed Ellipsoidal Dipoles 



Here is a mathematical puzzle of possible chemical interest. It is a relative of a popular statistical 
mechanics model, the dipolar liquid [27] . Consider an infinite set of hard ellipsoids of revolution (radius 
of revolution a/2, axial half-length 6/2. Let there be a permanent moment /i along the unique axis. 
What is the ground state structure and energy? This section only states and does not pretend to solve 
the problem, but makes a conjecture about the answer for the special case a = b of hard-sphere dipoles. 
This puzzle would be relevant to liquid phase self-assembly of ellipsoidal shaped molecules, if no other 
interactions were important besides hard-core repulsion and electrostatic interactions. 

One possible family of structures is arrays of "chains," where the term "chain" will be used to denote 
a stack of touching ellipsoids with the axis of revolution the same as the axis of the chain. Dipole- 
dipole interactions will favor ferro alignment along the chain. What transverse arrangement of chains 
is favored? Close packing will require at least an approximate triangular pattern of chain axes. In the 
"simple hexagonal" structure of perfect triangular symmetry and zero longitudinal offsets, first neighbor 
chains all wish to have antiparallcl moments. This is geometrically impossible, so the pattern is frustrated. 
Various anti-ferroelectric (e.g. alternating rows of up and down oriented chains) or ferrielectric (e.g. one 
out of three anti-aligned chains in a hexagonal pattern) give two out of three neighbors favorably anti- 
aligned and one out of three unfavorably aligned. 

The "offset" is defined as the distance along the chain axis separating centers of ellipsoids on two 
neighboring chains, measured in units of the period b which is also the axial length of the ellipsoid. With 
offset zero, touching ellipsoids have transverse separation r = a. With non-zero offset, the separation 
is r = ay/1 — (offset) 2 . Introducing offsets allows two good things to happen. (1) The chains can pack 
more densely because the widest parts of adjacent ellipsoids are not touching. (2) For offsets roughly 
between 1/4 and 1/2, parallel moments arc preferred. If all neighbors prefer to be parallel, there is no 
frustration problem. Of course, farther neighbors have other offsets, and if b/a > 1, farther interactions 
acquire more importance, reintroducing frustration. 

From the discussion of the previous section, two symmetrical possibilities can be seen: (1) a triangular 
array with offsets ±1/3, which gives the rhombohedral family to which the primitive cubic structures 
belong; or (2) a body-centered orthorhombic (bco) family to which the fee structure belongs. These 
structures are shown in Fig. 6. The bco structure minimizes the frustration problem by breaking triangular 
symmetry, giving the four closest neighbors the optimal offset 1/2 for ferroelectric order, and having two 
slightly more distant neighbors with zero offset which are unhappily ferroelectric. 

It is easy to see that bco is the preferred structure for dipolar hard spheres. Note that the density of 
ellipsoids (1/v where v, the volume per ellipsoid, is a 2 b/\/2 for bco and ^/l6/27a 2 6 for rhombohedral) is 
higher for bco at all axial ratios b/a. For cubic structures (these structures are fee and bee respectively 
when b = a) the CM law tells us that the dipole energy E = —(2n/3v)^ 2 depends only on the density. 
For other axial ratios, the CM law no longer works because the structure is not cubic. One might expect 
that bco would always win because of its higher density. However, numerical calculation shows that for 
axial ratios b/a> 1.717, the rhombohedral structure is better. This numerical test was done only in the 
range 0.5 < b/a < 2.5. For larger b/a, the on-chain parallel alignment of dipoles is no longer a big part 
of the energy, so there is no reason to think that chains would form. For smaller b/a, the chain-chain 
interaction is no longer a big part of the energy, so it is not likely that dipole interactions would play the 
dominant role in determining how chains arrange themselves. 

It seems reasonable to conjecture that the optimal solution for hard sphere dipoles is the bco=fcc 
structure with minimum volume. However, the complexity of dipolar interactions is such that it is not 
wise to be to confident of this solution. For ellipsoids that are close to spherical, the same answer probably 
applies, but for b/a = 1.717 where bco loses to rhombohedral, there is no good argument preventing a 
totally different structure from doing even better. 

9 Metallic Coatings 

The CMM focusses on a single aspect of the physics of electrical polarity in nanosystems, namely dipolc- 
dipole interactions of induced moments. In this final section the model is pushed a little further to 
model the effect of metallic screening from adjacent metallic regions. An interesting result was derived 
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Figure 6: Rhombohcdral (left) and body-centered orthorhombic (bco, right) packing of hard ellipsoids. 
The axis of revolution is shown vertical. Rhombohedral packing achieves a greater 2-d density of ellipsoids, 
but lesser 3-d density, than bco. 



by Maschhoff and Cowin [36] . They find that a neutral molecule with both a fixed point dipole moment 
^io and an isotropic polarizability a, when favorably aligned perpendicular to a metal surface at distance 
6/2, has energy (attraction to the image dipole) 

E = -» 2 /(b 3 - 2a). (20) 

The induced dipole enhances the total dipole and increases the attractive interaction energy between the 
dipole and the surface. At a critical polarizability a c — b 3 /2, the response is divergent. At this point, an 
induced dipole would spontaneously form even if the permanent moment fi were zero. 

It is evident that coupling a polarizable medium to metallic electrodes or shielding devices may 
promote spontaneous polarization. As a specific model, we now investigate an array of N fixed polarizable 
molecules (positions Ri) which are coated, or otherwise mingled, with N' fixed chargeable molecules 
(positions Re). The indices for polar sites are i, j and for charge sites are £, m. Many models of this type 
can be found in the literature, for example [7, 8, 9, 10, 11, 17]. The energy of chargeable sites alone would 
be the classical electrostatic energy 



i N ' i 



,ext 



2C i 2 e% Rem 

= ±<q\l(l/C)l + V}\q> + <q\^t> ■ (21) 

The term qf/2C is the energy cost of putting charge qe at the £'th site. Here \q > is an TV'-vector and 
1 and V are N' x N' matrices. The external potential <^ jCXt = 4> C xt{Re) is the external scalar potential, 
related to the external electric field by F cxt = — V<j> ext (r). The charges can adjust to minimize energy, 
subject to the constraint of charge conservation, qe =const, where the constant will normally be 
zero. Introducing a Lagrange multiplier A, we find that for a system of charge sites in constant external 
potential, 

[{l/C)l + V]\q>+\<j> ext >=A|0> (22) 

where < i\q >= qe is the charge at site i and < ^ 1 >= 1 defines |0 >. The first term on the left hand 
side of Eq.(22) is the vector \<\>- m t > of internal electrostatic potentials at sites. The charges adjust to 
make a constant potential <j> = </>; n t + <j> cx t- A non-zero A just uniformly shifts the potential in order to 
adjust the total charge. The system of charges should be neutral when not driven. This requires the 
matrix [(1/C)1 + V] to be positive. The Coulomb interaction V by itself is not positive - two sites with 
charge +Q and — Q have negative energy - so a capacitance C less than infinity is needed to keep charges 
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from spontaneously separating. In the resting state (</> e xt = 0) we want a neutral system (< 0\q >= 0) 
so A = 0. 

Now introduce electrostatic coupling between the charge sites and the polar sites, 

charges dipoles ^ _^ 

£int{{qe},{m}) = V] qi p 3 M ' 

= < q\K\n > (23) 

where the matrix K has dimension N' x 37V. For simplicity, choose <fi cxt — 0, and find the charge 
distribution \q > which minimizes the total energy, the sum of Eqs.(3,21,23), for an arbitrary fixed dipole 
distribution \fi >. The answer is 

1 ^ 1 



\i>=-[c + v ) K ^ > - (24) 

When this answer is inserted into the total energy formula, the result is an effective dipole-dipole inter- 
action r e ff, 

£tot = ^<M|[(Va)l-reff|M> (25) 

r eff = r + A't[(i/c)i + v]- 1 A: (26) 

where K* is the transpose of the rectangular matrix K . The charge system effectively enhances the 
dipole-dipole interaction, driving the system closer to a ferroelectric instability. 

As an illustration, consider a single polarizable molecule, with two chargeable sites, one to each side 
at a distance r. Assuming that charge responds instantly to polarization on the polarizable molecule, the 
energy to create a dipole is ^ 2 /2a ff, where 

J_ = I ±£—. (27) 

a eS a 1/2C - l/4r V ' 

The 1/r 4 in the numerator is the two powers of the charge dipole interaction K, and l/4r in the denom- 
inator is half the Coulomb attraction between the induced charges ±q. When there is an external field 
(along the trimer axis) the total dipole fi to t — H + 2gr has both a contribution from the polar molecule 
and from the metal, and the response is 



l/a-2/r 3 

1 H 

1 fi- - ±) 



F cxt (28) 



This diverges at the same critical polarizability that makes a c g diverge. 

It is easy to produce an alternative theory in which charges are treated by quantum mechanics, and 
charge interactions are solved in random phase approximation. The simplest such model is to let the 
charge sites Re be atoms which have only an s-orbital available near the Fermi level. The Hamiltonian is 

Hd = ^™4 c m + \ Vi™c\ c i c ™ c m ( 29 ) 
where c m destroys an electron from the orbital on site R m . The charge-dipole interaction is 

^el— dip ^ ^ -^m,ia/^zacJ n C m (30) 
iam 

Now the induced charge at site £ is related to the external potential at site m by the density response 
function x, 

Sqt = -^2 Xe m SV extim . (31) 
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and the energy of the induced charge distribution (analog of Eq.(21) is 

£c*({qi}) = \<q\x- 1 \q> (32) 

If we ignore the Coulomb interaction (second term in Eq.(29)), then the electrons are governed by the 
non- interacting Hamiltonian H® v This problem has the solutions H^n >= e n \n >, and the susceptibility 
can be computed from 

occ empty . , , . . . 

o _ y^Y^ < £\n >< n'\m >< m\n X n\£ > 
Xe m - 2 2_^ ; — ~ ( dd > 

n n> n " 

However, it is crucial to account for the Coulomb interaction between electrons, which means solving the 
electron many-body problem. This cannot be done exactly. The random-phase approximation (RPA, 
also known as the time-dependent Hartree approximation) provides good insight and not a bad numerical 
answer. The relevant equation is 



Then Eq.(32) is replaced by 



Xl^ext >= Xod^ext > + 1 (frVkartree >) (34) 

£*({<//}) « \ < llXo 1 + V\q > (35) 

Note that the structure of the classical model Eq.(21) is closely followed, with the non-interacting sus- 
ceptibility interpreted as a non-local capacitance. The result is an effective dipole-dipole coupling which 
replaces the classical Eq.(26), 

r cff = r + ^t[i/ Xo + T/]- 1 ^ (36) 

As an illustration, consider the same problem of a single polarizable molecule, with two chargeable 
sites, one to each side at a distance r. This time the sites have s-orbitals coupled by hopping matrix 
element t and Coulomb coupling q 2 /2r. It is also sensible to include an on-site Hubbard repulsion U 
which discourages deviation from charge neutrality. This system responds to external fields with the 
effective polarizability 

1 1 1/r 4 
_!_ = til (37) 

a eS a t/e 2 + U/2e 2 - l/4r" V 1 

Halas et al. [37, 38] have discovered how to coat semiconducting nanoparticles with metal shells. It 
would be interesting to attempt this with nanoparticles of large polarizability, in order to engineer new 
superparaelectric or nanoferroelectric systems. 

Note added in proof, Nov. 15, 2003 

Concerning the possiblity of fcrroelectricity in fullcrenes via the CM polarization catastrophe, R. L. 
Whetten had suggested to me endohedral fullerenes as a route to enhanced polarizability. It turns out 
the Clougherty [39, 40] has speculated about ferroelectricity in such systems, not from a CM point of 
view, but more correctly, from a Jahn- Teller point of view. 

N. W. Ashcroft has reminded me that the CM polarization catastrophe was formulated by Goldham- 
mer [41] and Hcrzfcld [42] as a criterion for metallization under pressure. As the macroscopic dielectric 
function diverges, we now understand that the frequency of some "ferroelectric soft mode" goes to zero. 
If this is not a lattice mode, but an electronic mode (plasmon), then the new phase as understood by 
Goldhammer and Herzfeld, should be a metal, not a ferroelectric. This criterion proved useful for guiding 
high pressure experiments [43], particularly for the metallization of Xe [44, 45], where there are no lattice 
modes internal to the "molecule" of Xc. 

N. W. Ashcroft has also reminded me of an interesting extension of dipole ordering energy theory, 
namely the construction of harmonic polarization wave states, first done by Lundqvist and Sjolander 
[46]. Atwal and Ashcroft are applying this idea to superconducting interactions [47], following earlier 
work [48]. 

Finally G. Stcll has introduced me to some of the literature on computer simulations of dipolar systems. 
Some of these papers [49, 50, 51, 52, 53, 54, 55, 56, 57] are particularly relevant to the discussions in the 
present paper. 
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